Exact results for nonlinear ac-transport through a resonant level model 
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We obtain exact results for the transport through a resonant level model (noninteracting Anderson 
impurity model) for rectangular voltage bias as a function of time. We study both the transient 
behavior after switching on the tunneling at time t = and the ensuing steady state behavior. 
Explicit expressions are obtained for the ac-current in the linear response regime and beyond for 
large voltage bias. Among other effects, we observe current ringing and PAT (photon assisted 
tunneling) oscillations. 

PACS numbers: 



I. INTRODUCTION 

The recent advances in nanotechnology created a lot 
of interest in transport through correlated quantum im- 
purities. While the linear response regime essentially 
probes the ground state properties of the system, trans- 
port beyond the linear response regime explores genuine 
non-equilibrium quantum many-body phenomena. How- 
ever, theoretical calculations beyond the linear response 
regime are challenging since the steady state cannot be 
constructed via a variational principle like equilibrium 
states. Even for dc-bias only recently exact numerical 
methods have been developed that permit such investiga- 
tions for interacting systems, notably the time-dependent 
numerical renormalization group [1 , Monte Carlo meth- 
ods [2, 3 , and the time-dependent density matrix renor- 
malization group [H |5]. Some of the analytical meth- 
ods that have been applied successfully are perturbative 
Keldysh calculations [6 , extensions of the renormaliza- 
tion group [71 [8] , flow equations , and generalizations 
of NCA (non-crossing approximation) to non-equilibrium 
[10l[TT] . A comparative review of theoretical methods can 
be found in Ref. [12 

For ac-bias beyond the linear response regime still 
much less is known since, e.g., the numerical methods 
cannot easily be generalized to time-dependent bias. In- 
teresting ac-phenomena are for example the photon as- 
sisted tunneling effect (PAT) [13] that has been observed 
in experiments , or the " current ringing" after a step- 
like bias puis p31. Non-equilibrium Green's function 
methods can be employed p!5HT7] when the correlation ef- 
fects are not too strong. In the strongly correlated regime 
of the Kondo model the non-crossing approximation was 
found to be reliable JT] IT8H2Q] . At a specific point of 
the two-lead Kondo model it can be solved exactly [21 , 
which permits exact results for the current in the steady 
state [22j, after a rectangular pulse 22 under sinu- 
soidal bias [2T1. Unfortunately, this special point is not 
generic for a Kondo impurity that can be derived from 
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an underlying Anderson impurity model, which is exper- 
imentally the most relevant situation. 

In this paper we study the response of a resonant level 
model (noninteracting Anderson impurity model) under 
rectangular ac voltage bias after switching on the tunnel- 
ing at time t = 0. We derive exact analytical results for 
the transient and steady-state current by diagonalizing 
the Hamiltonian. This exact solution contains both dc- 
and ac-bias in and beyond the linear response regime. 
While dc-results and ac-results with sinusoidal bias have 
been obtained previously in the literature [15] , rectangu- 
lar ac-driving beyond the linear regime seems not to have 
been studied before. Besides being experimentally rele- 
vant, our results are also helpful for exploring the various 
crossovers in this important model and serve as an exact 
benchmark for future work. 



II. MODEL AND DIAGONALIZATION 

The resonant level model coupled to two leads is de- 
fined by the following Hamiltonian 

H = ^efc4^cfc« + ^-^(4^d+/i.c.), 

ka koi 

where a = L^R denotes the leads. The spin index can 
be omitted since the model is non-interacting and we 
work with spinless fermions. All energies are measured 
with respect to the single-particle energy of the impurity 
orbital (e^^ = 0). We take a wide band limit with a linear 
dispersion relation, Ck = kr]^ where r] denotes the level 
spacing and k an integer number. The hybridization is 
defined by F = piig^ where p = l/rf. The impurity orbital 
spectral function in equilibrium is then given by 

p^i^) = ^^g^^ (1) 

Our strategy to obtain exact results is to first diagonal- 
ize the discretized Hamiltonian and to then take the ther- 
modynamic limit 77 ^ 0. We introduce the hybridized 
basis Cs = X]/c e -efc ^sCk-\- + Bsd. It is then straightfor- 
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ward to diagonalize the Hamiltonian 



C„Ca 



(2) 



where Ck± = -^{ckL i Ckn)- The inverse transformation 
is d = BsCs and c/e+ = j^^BsCs. The eigenval- 



ues are determined as solutions of the equation 
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In the thermodynamic hmit 
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From the diagonahzation one also derives the following 
set of equations 
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which will be important for calculating various summa- 
tions below. 

An ac voltage bias leads to time-dependent potentials 
Ua{t) in the leads and the Hamiltonian takes the form 

^ = E(^fc " «"(t))cLcfea + ~/^'^^kad + h.c). (5) 



ka 



We suppose that initially (at time t < 0) the left and right 
lead chemical potential are the same, jul = fJ^R = the 
hybridization is switched off and that there is no elec- 
tron in the dot, Ud = 0. At time t = the hybridization 
is switched on and a rectangular voltage bias with pe- 
riod 2T (see Fig. 2) is applied: UR{t) = -UL{t) = V/2 
for 2NT < t < (27V + 1)T and UR{t) = -UL{t) = -V/2 
for {2N + 1)T < t < 2{N + 1)T.[^ fi therefore gives the 
energy difference of the impurity orbital to the " average" 
Fermi energy of the leads for time t > (Fig. 2). 
The current operator is defined as 



k 



V2 



where denotes the total number of electrons in lead a 

-, def def 

and sl = I, sr = —1. 

In the first half period 2NT < t < {2N + 1)T the 
Hamiltonian is 



Ha = ^{^k^^)clLCkL^^{ 

k k 
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0<t<T 



T<t<2T 





FIG. 1: A schematic diagram of our model: A step- like volt- 
age bias is applied to the two leads coupled to the quantum 
dot. 



Because the dispersion relation is linear and k runs 
from — oo to oo (wide band limit), we can simply re- 
label the fermion operators, Cka = q,- The po- 
tentials in the leads are eliminated by this transforma- 
tion and the Hamiltonian can be diagonalized as be- 
fore: Ha = Y^s^sdlds ^ Y.k^kCi\_OLk-, where = 



E 



-ak+ + Bsd and ak± = -^{cj. 



k TJ^^k^ -r J^s^ ciii^ a./e± - y^V'-/e-^,L ^ ^k+^^.RJ- 

Similarly, in the second half period (2A^ + 1)T < t < 
2{N + 1)T the Hamiltonian is diagonalized as = 



and bk± = 4^(c, 



-Y.k ^kh\_hk-, where hs = Yl 



In the Heisenberg picture the current operator at time 
t = 2NT-\-r^ r G [0, T] (first half period) can be expressed 



k es-ek 



-B,d 



and in the second half period {t = {2N + 1)T + r, r G 
[0,T]) 



{d^Cka-clj). (6) e.g., Es^y 



xe ^^^.^'(e-^^^'^e-^^^^)^\ (9) 



To find Ia{t) we first calculate the time evolution of the 
single fermion operator d^ and c^^ under Ha or by 
expressing d^ and c^^ in the hybridized basis, next apply- 
ing the diagonal time evolution and finally transforming 
back to the original basis. The calculation is straightfor- 
ward but one needs to pay attention when encountering 
summations with respect to the eigenenergies e^. In the 
thermodynamic limit the summation can be transformed 
into an integral when there is no pole in the integrand, 
/ despBle~^^^^ = e~^^ If there are 



poles in the integrand we first calculate the time deriva- 
tive to get rid of the pole terms. Key formulas are 



-icst 
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By using these two formulas we get ^^^^^ ^(a,b) ^ , = = + V/2 

e^H^a,b)T^^ -iH^a,b)T ^ g-rT^t _^ ^ ^ Vl/^'''^^c^^l2) aiid e^^ = e^^ = e/, - y/2. Employing this formula twice 

gives the evolution over a full period: 

(13) 

I 



where 



I 

We perform the summation over k by transforming it Applying the above formula recursively times yields 
into an integral and then employing the residue theorem. 



iHaTiHbT\Nj {^-iHtT-iHaT\N 
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where ao = /3o = 70 = and the recursion relations are 



Q^iv+l(A:^ A:) 
^iv+l(A:^A:) 
lN+i{k) 

It is easy to find 
In 
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In the first half period the current evaluates to 
h 



-2Im(e, 



0-) g-2NiekT-iel 



where ^i^j = ^L7iv(A:)e-r- ^ g2iVie,T|^a^(^)^ ^j^^ 
Fermi-Dirac distribution function. In the sequel we will 
always specialize to the zero temperature case {uk = 1 
for A: < 0, n/c = for A: > 0). In the second half period 
the current evaluates to 

(2) |2 



-2Im(ei'je-^^^^^^-^^^«^-^^'^«^)) (20) 
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where ^fj = 



e"^^fc«^iy^^(r)). To simplify notation in lengthy expres- 
sions we will frequently employ F as the unit of energy 
and current, and l/F as the unit of time. In the final re- 
sults we always reintroduce all dimensionful parameters. 



2NiekT 



-Ft 



III. BUILDUP OF THE STEADY STATE 

There is a transient time regime after the coupling of 
the dot to the leads is switched on at time t = before a 
steady state has built up. Initially, the left lead current 
is opposite to the right one and the initially empty dot 
is being charged. We will see that these transient effects 
decay exponentially (proportional to e~^^) to the steady 
state. 

Let us explicitly look at the two limits of period T 
oo (dc bias) and T ^ (very fast driving). For T ^ oo 
one finds from Eq. (19) 
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FIG. 2: Time-dependent current for different switching pe- 
riods T of the ac voltage bias (top: infinitely fast driving, 
middle: intermediate fast driving, bottom: dc case). Zero 
temperature and ac voltage bias V — V in all graphs. The 
full lines denote the left lead current, the dashed lines the 
right lead current. The hybridization is switched on at time 
t = 0. Notice the discontinuous onset of the current at t = 0, 
which is due to the wide band limit for the conduction band 
(a detailed discussion can be found in Ref. [26 ). 
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The steady limit {t oo) is 
eT f , , , eV 
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which of course coincides with the well-known result for 
the stationary dc-current [15 , e.g. for zero temperature 



2eF feV 
— — arctan — — 
h V 2F 



(23) 



In the fast driving limit T ^ we keep t = 2NT 
invariant and let ^ oo. According to the Trotter 
formula, the evolution then becomes equivalent to zero 
voltage bias [25 , limT^o(e^^^^e^^^^)^ = e^(^-+^^)™. 
We find 



-^a (^) — 



2eTe- 
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In Fig. 2 we show the transient currents in the left and 
right lead for different periods T when /i = 0. The cur- 
rent oscillations are suppressed when the frequency goes 
to infinity. The /(t) -curves gradually change from the dc 
limit to the high frequency limit described by Eq. (24) 



when the period T decreases. In the fast driving limit 
the left and right currents are opposite to each other and 
both decay to zero with increasing time. 



IV. STEADY STATE BEHAVIOR 

When the time is much larger than l/F, the current 
reaches its steady state behavior. By taking A/" — > oo we 



find this steady state limit given by 

eT f ~ ~ 

Ia{r) = Sa-^ j deknk{\£,kLf + - 2lm^ka), (25) 

where < r < T. In the first half period we have 

1 ScKV(e^*^^~*^fc«^~^ — e^^ko.iT-^)-T—^^ 



2ieT 



(26) 



and in the second half period 

1 s„T/(e«L(r-r)-T-T _ g2ieT-ieLr-r-, 



(27) 



From Eqs. (26) and (27) one immediately verifies that 
the steady state current satisfies /^(r) = — /^(r + T) as 
expected intuitively, where a denotes the opposite lead. 



A. Linear response regime 

In the linear response regime of small voltage bias a 
sinusoidal signal drives a sinusoidal current with the same 
frequency, and signals with different frequencies can be 
superimposed linearly. Therefore we can factorize the 
rectangular signal into a series of sinusoidal components 
and find the frequency-dependent complex admittance of 
the system. 
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FIG. 3: The linear admittance of a resonant level model for 
various level positions /x (energy of the dot level with respect 
to the Fermi energy in the leads) at zero temperature. The 
top graph shows the absolute value of the admittance, the 
bottom one its phase. 



In the linear response regime the left lead current is 
equal to the right lead one and can be expressed as 



lim 



V 



den{e)T{e), 



where 



T(e) = 



2er3 
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We Fourier transform both the ac-voltage signal and 
the current. We define I{uJn) = J^^ dte'^^'^^I{t) = 
2 dre^^''^ I{t)^ where we use the property /(r + T) = 

-/(r), and V{uJn) = J^^ dte'^^-^Vit). Here uJn = ^ 
and n is an odd number. The voltage bias is —V 
for < t < T and V for T < t < 2T, leading to 
V{uJn) = By adjusting T the frequency uj^ can be 
an arbitrary real number, and the linear response admit- 
tance G{uj) = I{uj)/V{uj) at zero temperature is given 

by 



Giu;) = 



I arccot- 



arccot- 



h 



2uj/T 



2^2 



-F2)((/i-c^)2+F2) 



(30) 
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where /i denotes the position of the dot level with re- 
spect to the average Fermi energy of the leads, see Fig. 1. 
Eq. (30) agrees with previous ac-calculations in the lin- 
ear response regime, see Ref. [27^. Fig. 3 depicts G{uj) 
for different level positions /i. The admittance goes to 
zero for fast driving, a; ^ oo. For a; ^ one recovers 
the well-known dc-conductance G = ^ ij,'^-\-r'^ • -^^^ asym- 
metric dot positions the resonance peak is around uj = ji^ 
showing the PAT (photon assisted tunneling) effect [M]: 
When the frequency of the ac-signal is equal to the en- 
ergy difference of the dot level from the Fermi energy in 
the leads, electrons in the leads can absorb a photon and 
jump into the dot. Notice from Fig. 3 that the symmetric 
dot always acts like an inductor as already explained in 
Ref. [27]. For asymmetric dots there is a crossover from 
capacitive to inductive behavior around uj = ja \27\ . 



B. Beyond the linear response regime 

For a voltage bias beyond the linear response regime it 
is impossible to calculate G{uj) by performing a Fourier 
transformation since the different frequency components 
interact with each other nonlinearly. Therefore we now 
depict the behavior of the current I{t) as a function of 
time t during one full period in the steady state situ- 



ation. Due to the nonlinearities we need to discuss this 
separately for different driving periods T. We will always 
take zero temperature in the sequel, the generalization to 
nonzero temperature is straightforward. 

We first look at fast driving, T <C F~^. For the sym- 
metric situation the I — t curve becomes triangled: The 
current decreases from maximum to minimum in the first 
half period, and then increases from minimum to maxi- 
mum in the second half period, see Fig. 4. In the oppo- 
site slow driving limit T :^r~^, the I — t curve becomes 
rectangled. The saturated current in each half period 
is simply given by the corresponding steady dc-current 
(23). For intermediate driving speed, T ~ F"^, we ob- 



serve ringing oscillations [15] of the current with period 
Att/V (see Fig. 5). 

For asymmetric dot positions /i 7^ the current also 
has characteristics of PAT and ringing, which are, how- 
ever, not easily visible in a plot like Fig. 6. Clear sig- 
natures can be found in the the differential conductance 
with respect to the gate voltage, which we denote as gate 
differential conductance G^ate distinguish it from the 
usual definition of differential conductance with respect 
to the voltage bias between the leads. We define 



def dia (r) 



and the current can then be expressed as /^(r) 



(31) 
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FIG. 4: The steady state current in one period for fast driv- 
ing (here T = O.l/F and zero temperature) in a symmetric 
resonant level model (/i = 0). The left figure depicts the cur- 
rent / (units eV/h) in the nonlinear regime, and the right 
figure shows I /V (units /h) for smaller voltage bias (linear 
response regime). Because the driving period is shorter than 
the time required to establish stationarity in one period, the 
time-dependent current looks triangular. 
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FIG. 6: The steady state current in one period for inter- 
mediate driving (here T = 1/F and zero temperature) in an 
asymmetric resonant level model (/x = 5F). The left figure de- 
picts the current / (units eV/h) in the nonlinear regime, and 
the right figure shows I/V (units e'^/h) for smaller voltage 
bias (linear response regime). The crossover from capacitive 
to inductive response (compare Fig. 3) leads to a complicated 
behavior of the current in the first half period. 



Current I (eF/h) 



1^60 



•^40 



o20 




W (e^/h) 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



T 

Time t 



2T 




Time t 



Time t 



FIG. 5: The steady state current in one period for inter- 
mediate driving (here T = 1/F and zero temperature) in a 
symmetric resonant level model (/i = 0). The left figure de- 
picts the current / (units eV/h) in the nonlinear regime, and 
the right figure shows I/V (units e^/h) for smaller voltage 
bias (linear response regime). The oscillations of the current 
with period Atv/V ("current ringing" J^) are clearly visible 
for large bias. 



Figs. 7 and 8 shows G^^*® in the first half period 
(^Qgate second half period follows via G^fnd^(e) = 

Gfg^^(— e)). In the linear response regime we find a pair 
of bright PAT lines at e = ±7r/T (see Fig. 7). In the 
regime far from equilibrium, high order PAT lines at 
e = nn/T {\n\ > 2) can be observed (see Fig. 8), indicat- 
ing multiple photon assisted tunneling processes. These 
PAT lines combine and are replaced by a pair of bright 



FIG. 7: The gate differential conductance Gj^^^{e,T) (units 
/h) in the linear response regime [V — 0.2F and zero 
temperature) for period T — 0.2/F in the left figure and 
T — 0.6/F in the right figure. The pair of bright lines sym- 
metric to e = are the PAT lines at e = ±7r/T. 



resonance lines at e = ±y/2 when the period increases. 
This demonstrates that ac transport for high frequencies 
is dominated by photon assisted tunneling, and by reso- 
nance tunneling for low frequencies. 

V. CONCLUSIONS 

We have investigated a resonant level model driven by 
rectangular ac-bias in and beyond the linear response 
regime. Even this simple model shows surprisingly rich 
behavior in its transport properties. One can observe 



7 




specific nonequilibrium effects like the buildup of the 
steady state, current ringing and photon assisted tunnel- 
ing, and the crossover to the well-studied limiting cases 
of dc-bias and linear response regime. The results are 
exact and based on an explicit diagonalization of the 
Hamiltonian in the first and second half period of the 
rectangular voltage bias driving. Within the flow equa- 
tion framework, this approach can easily be generalized 
to an interacting quantum impurity model exposed to 
ac-driving beyond the linear regime. Much less is known 
about such systems, which provides another motivation 
for this work and will be studied in a subsequent publi- 
cation. 



FIG. 8: The gate differential conductance Gf 



(units 



e^//i) for large voltage bias (V = 20r) and zero temperature. 
The top left figure shows fast driving (T = 0.2/r), T = 0.5/r 
in the top right figure, intermediate driving (T = 1/r) bottom 
left and slow driving (T = 5/r) bottom right. The y-axis 
denotes the energy ranging from — 20r to 20r. For fast driving 
(T < 0.8/r) the higher-order PAT lines are clearly visible. For 
slower driving (T > O.S/F) the PAT lines away from e = ±V/2 
disappear with increasing T. 
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